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ABSTRACT 

We study the UV and optical properties of 38 massive globular clusters (GCs) 



m 
o 

o 

O in the Virgo elliptical, M87, imaged using the ST IS and WFPC2 instruments 

O ' onboard the Hubble Space Telescope. The majority of these GCs appear ex- 

| tremely bright in the far-ultraviolet (FUV) - roughly a magnitude brighter 

than their Galactic counterparts with similar metallicities. The observed FUV 
flux is several times larger than predictions of canonical old stellar population 

X! ■ 

' models. These canonical models, which assume a fiducial helium enrichment 

parameter, AY/ AZ — 2, are able to reproduce the observed FUV fluxes only 
if ages ^3-5 Gyr larger than the "WMAP age" of the Universe are invoked, 
although the same models fit the UV photometry of Galactic and M31 GCs 
for ages less than the "WMAP age". A similar discrepancy (~ 3 Gyr) is found 
between the mass-weighted and f/V^luminosity weighted ages of the massive 
Galactic GC u> Cen, whose colour-magnitude diagram (including peculiar fea- 
tures on its well-populated horizontal branch) can be accurately reproduced 
by invoking a small super-He-rich (AY/ AZ > 90) stellar component. By com- 
parison to uj Cen, we propose that the majority of M87 GCs in our sample 
contain strong signatures of similarly minor super-He-rich sub-components. 
This hypothesis is supported by simulations which suggest that, based on 
the UV detection limit of this survey, the number of GCs detected is several 
times of the prediction from canonical models. Although we cannot prove or 
disprove the extreme helium scenario at the moment, we show that the same 
phenomenon that causes the extended horizontal branch of u Cen explains 
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the UV brightness of our sample. If this is indeed due to the extreme helium, 
this study would be the first to find its signatures in extragalactic objects. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - 
galaxies: formation - galaxies: fundamental parameters - globular clusters: 
general 



1 INTRODUCTION 



Recent observations of the stellar population of the largest Galactic globular cluster (GC), 
L) Cen, have revealed a double main sequence (MS), with a disjoint sub-population of bluer 
and fainter MS stars which is sep arated from a majority population of redder and brighter 
MS stars (e.g. iBedin et al.l 12004 ). In addition, ui Cen harbours a substantial population 



of extended horizontal branch (EHB) stars, which are hotter and distinct from the primary 



population of normal hor izontal branch (HB) stars (jWhitney et al 



1994 



Ferraro et al 



Sollima et al. 



2005 



2004 : 



20061 ). Theoretical analysis has indicated that the double main se quence 



could be a result of a large spread in helium (He) abundances within the cluster (INorris 



20041 ) - this is supported b y spectroscopic measurements o f metal abundances of MS stars in 



uj Cen ( Piotto et al. 



20o3). More recently. 



Lee et al. 



(120051 ) have shown that a large variation 



in the He abundance in uj Cen can not only reproduce the double MS, but also accounts 
naturally for the 'anomalous' population of hot extended HB stars in this clustej^]. Their 
best simulation of uj Cen requires the He enrichment (AY/AZ) parameter, for a minor- 
ity populat ion of stars, to be in excess of 90 - for comparison AY/AZ ~ 2 on galactic 



Peimbert et al. 



scales (e.g., 
ficult to reproduce theoretically (e.g., 



20011 ). Alt hough such hig 



i values of helium enrichment are dif- 



Bekki fc Norris 



2006 



Choi fc Yi 



20071 : 



Karakas et al 



20061 ). similar super- He-rich sub-populations are also capable of reproducing the peculiar 



HB morphology of another Galactic GC, N GC 2808 (ID'Antona fc Calo 



2004 



Bedin et al 



2000 



Lee et al. 



2005 



D'AntonaetaL 



20051 ) . A similar phenomenon may be responsible for 



the unusually-blue HB morphology of NGC 6338, a nd NGC 6441 for their metallicities (e.g., 



Sweigart fc Catela: 



1998 



Caloi fc D'AntonalbOQfih . 



Helium- rich populations are thought to have a significantly higher number of hot, evolved 
HB stars and thus show enhanced UV fluxes (and hence bluer integrated UV colours) 

* Send offprint request to yi@yonsei.ac.kr 

Note that the evidence for He enrichment in u Cen does not come from direct observations of He spectral lines but is inferred 

from the HB morphology and direct determinations of metallicity. 
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than helium -poor counterparts of the same age (e.g. 



1998 



Dorman et al. 



1995 



Yi et al. 



19971 . 



19991 ) . Using models with standard He enrichment (AY/AZ ~ 2) to derive ages 



from such integrated UV colours would therefore lead to those ages being overestimated, 
since populations with standard He enrichment take longer to output the same level of UV 
flux as those with a super-He- rich sub-component. Reversing this argument implies that 
the detection of GCs that are too blue to be produced by old populations with standard He 
enrichment, might indicate the presence of anomalously high He concentrations in the stellar 
populations within these clusters. In this paper we present a plausibility argument for the 
presence of just such a population of GCs in the elliptical galaxy M87. The motivation fo r 



Sohn et al.l 



2006|). 



this work stems from the recent UV study of M87 globular clusters (see 
We find that (a) the FUV detection rate of M87 globular clusters in this study is significantly 
larger than can be expected from simulations of a standard old stellar population that has a 
He enrichment value of ~ 2 and that (b) the majority of globular clusters which are detected 
have derived (FUV -weighted) ages which are far in excess of the presently-accepted age of 
the Universe. In this study, we adopt the most rec ent estimate from the WMAP data as the 
"physical" age of the Universe (ISpergel et al.ll2006l ). Since the stellar models used to perform 
the GC age estimation are well-calibrated to Galactic GCs, and produce physical ages (< 13 
Gyr) for old GCs in M31, we argue that both pieces of evidence are consistent with the 
existence of super-He-rich populations in the GC system of M87, and that indeed it may 
be the FUV -enhanced 'super-helium-enriched' tail of the GC distribution which has been 
sampled by this study, given that the FUV detection rate is far higher than that expected 
for a normal population of HB stars. 



2 PARAMETER ESTIMATION: AGE AND METALLICITY OF M87 GCS 



We be gin by combining HST '/ ST IS broad-band FUV and NUV photometr y from 



Sohn et al. 



Jordan et al 



fj2002h i n the F336w, 



Kundu et al. 



(120061 ) with HST/WFPC2 optical photometry from 

F410m, F467m and F547m filters and photometry in the V and I filters from 
( 119991 ). to simultaneously estimate the ages and metallicities of 38 GCs in M87. GC star 
formation histories (SFHs) can be adequately parametrised by (dustless) simple stellar pop- 
ulations (SSPs) of a given age and metallicity alone, and our parameter estimati on em ploys 



a finely-interpolated grid of SSPs using the Yi-Yoon stellar models (see e.g. 
parameter estimation proceeds as follows: 



Yi 



2003|). The 
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For a vector X denoting parameters in the model, and a vector D denoting the measured 
observables (data), 

prob(X|D) oc prob(D|X) x prob(X), (1) 

where prob(X|D) is the probability of the model given the data (which is the object of 
interest), prob(D|X) is the probability of the data given the model and prob(X) is the prior 
probability distribution of the model parameters 

If we assume a uniform prior in our model parameters i.e. 

prob(X) = constant, (2) 
we have 

prob(X|£>) oc prob(£>|X). (3) 

Assuming Gaussian errors on the measured quantities gives 

probplX) oc exp(- X 2 /2), (4) 

where exp(— x 2 /2) is the likelihood function and x 2 is the sum of the normalised residuals 
between the model-predicted observables M k and the observed values D k 



k=i 



and <7fc is the error in the residual M k — D k . prob(X|D) is a joint PDF, dependent on 
all the model parameters. To isolate the effect of a single parameter Xi in, for example, a 
two parameter model (prob(X|L>) = prob(Xi, X 2 \D)) we can integrate out the effect of X 2 
to obtain the marginalised PDF for X x : 

roo 

probOTJD) = / prob(X 1 ,X 2 \D)dX 2 . (6) 
Jo 

In our analysis, the parameters X in the model are the age (t) and metallicity (Z) 
of the SSP in question, and the observables are colours constructed from UV and optical 
photometry. We choose a uniform prior in t, in the range 1 to 18 Gyr. For the metallicity 
Z, we use a uniform prior in the range —2.3 to 0.7 dex, which are essentially the metallicity 
limits of the stellar models. For each GC analysed, we marginalise the age and metallicity 
parameters, take the peak (x ) of the marginal PDF P(x) as the best estimate and define 
one-sigma limits x + and x_ as follows: 

/ P(x)dx = 0.34 / P(x)dx (7) 

Jxq JXq 
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Figure 1. Age and metallicity estimates for M87 GCs using UV and optical photometry (8 bands in total). 

/ P{x)dx = 0.34 / P(x)dx (8) 

We define the uncertainty in the estimate of a parameter as + i.e. the total extent 
of the one-sigma limits in parameter space. 



In Figure Q] we show the estimates of age and metallicity for each GC, derived using 
this procedure from a set of stellar models which has the fiducial value of He enrichment 
(AY/AZ = 2). We do not find a clear hint for a substantial di fference in age betw een the 



Jordan et al. 



(j2002|). 



metal-poor and metal-rich GCs, in agreement with the result of 

Figure [2] shows the marginalised probability density function (PDF) for age and metallic- 
ity for all M87 GCs used in this study. The metallicity PDF appears multi-modal. Applying 
the KMM test for bimodality yields two peaks: a metal-poor peak at —1.33 dex a nd a m etal- 
rich peak at —0.25 dex. These values agree r easonably wel l _ with I Jordan et al.l (120021 ). who 



report values of —1.58 dex and —0.30 dex and 



Kunduetal 



f|1999h . who find values of -1.41 



dex and —0.23 dex. The KMM test accepts bimodality at the 98 percent confidence level 
and trimodality at > 99.9 percent confidence. However, we note that the trimodality, in par- 
ticular, is weak and that the small size of the sample, coupled with its UV -selected nature 
makes it unreliable to extrapolate this metallicity PDF to the entire M87 GC population! 
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Figure 2. Marginalised probability density functions for age and metallicity for all M87 GCs used in this study. 

Since the metallicity estimation in this study is derived from 8 photometric filters instead 
of one single colour index, it is plausible that we are better able to resolve the metallicity of 
the population and, in particular, the bimodality of the metal-poor peak, which may not be 
obvious when using a single optical colour such as V — I. We cannot, however, estimate the 
proportion of each subpopulation from this data, because the GCs represent a FUV selected 
sample and are therefore unrepresentative of its parent population. 



3 THE DERIVED AGES OF M87 GLOBULAR CLUSTERS 

The ages derived in the last section, using the fiducial stellar models, range 14-18 Gyr, 
unphysically large given current estimates of the age of the universe. The technical reason 
for this apparent 'overestimate' of the age is that the M87 photometry largely lies outside 
the model colour-colour grids which cover the age of the universe (t ~ 14 Gyr). We illustrate 
this by showing a model (FUV — V) vs. (V — I) grid in Figure El While the Milky Way GC 
data (open triangles) are reproduced by the grid, most of the M87 data (filled circles) lie 
outside the age ran ge 1-14 Gyr for the generally-accepted metallicity range, [m/H]< +0.2 



(ICohen et al.lll998h . 
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Figure 3. Model (FUV — V) vs. (V — I) grid for a range of metallicities and ages, generated from stellar models with the 
fiducial value of He enrichment (AY/AZ = 2). The lowest age plotted is 1 Gyr and the largest age plotted is 15 Gyr. Ages 
1,5,10 and 15 Gyr are shown using filled circles of increasing sizes. The GC data of M87 (filled circles) and Milky Way (open 
triangles) with errors are ovcrplotted. It is apparent that the M87 photometry lies outside the age range 1-14 Gyr for all 
metallicities. 



It is important to note that the ages derived for the M87 GCs in the previous section 
are essentially FUV luminosity-weighted ages. For comparison, we now seek an estimate of 
the age of the M87 GC population which does not rely on their FUV flux. We note that 



Cohen et al. 



using the 



Worthey 



1998) have previously derived ages of M87 GCs (as a function of metallicity) 
(119941 ) models. We repeat their analysis, using their measurements of the 



Lic k indices of M87 GCs (see their Table 1 for the specific indices used), but this time using 



the 



Bruzual fc Chariot 



(120031 . BC2003 hereafter) stellar models, which are more up to date. 
Using the BC2003 models, our Lick index-based age estimates on Cohen et al.'s M 87 GCs 

(119981 ). Figure 



Cohen et al. 



are roughly 10-11 Gyr, which is somewhat lower than those of 
H] shows the same information as Figure [I], with t he ages determined from the Lick indices 



overplotted using open triangles. 



Lee et al. 



(120001 ) have previously noted that a population 
with a hot HB would appear younger in Lick indices by ~ 2-3 Gyr due to a significant 
contribution from the HB stars to the .B-band flux. Even considering this effect, the Lick 
index ages are still 2-4 G yr younger th a n the UV-dehved ages. 



We note here that the 



Cohen et al. 



(119981 ) sample does not overlap with the 



Sohn et al. 
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Figure 4. Same as Figure[T]but with the ages of M87 GCs (studied bv lCohen et al.l jl99Sl) ) determined from Li ck absorption 
line in dices and their uncertainties (open triangles with error bars and grey shaded region) added. See Table 1 of I Cohen et al.l 
dl998l) for the Lick indices used in this age determination. 



(120061 ) sample of GCs. Although we would ideally like to determine the 'spectroscopic' ages 
of the Sohn et al. GCs, no absorption line data exists for this sample and therefore we are 
not able to extract ages and metallicities for th ese objects indepen dent of their UV-optical 



colours. However, if we now assume that (a) the lCohen et al.l (119981 ) sample is representative 
of the parent M87 GC population and (b) the age derivation from Lick indices is closer to 
the true age of the M87 GC population than that derived from the UV photometry, then we 
find a ~ 2 — 3 Gyr discrepancy between the the FUV luminosity-weighted mean age of M87 
GCs in this sample and the 'true' age of the parent population. The discrepancy is driven 
by an enhancement of the FUV flux in these clusters beyond that expected from normal 
HB stars. 



4 FUV ENHANCEMENT IN M87 GCS: SUPER-HE-RICH STARS? 



Lee et al. 



(120051 ) have recently suggested that a small super-He-rich stellar component repro- 
duces both the peculiar features of the main sequence and the anomalous EHB population 
(which is hotter and emits significantly more FUV flux than normal HB stars) in u Cen. 
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Figure 5. Comparison between age and metallicity derivations using the fiducial (AY/AZ = 2, open red squares) and helium- 
enhanced (AY/AZ = 3, black crosses) stellar models. 



In this section, we explore the plausibility of a similar mechanism being responsible for the 
enhanced FUV flux seen in M87 GCs. 

In order to explore the effect of enhanced He enrichment on the age derivation of the M87 
GCs, we repeat the parameter estimation of Section 2 but replace the fiducial stellar model 
used previously (AY/AZ = 2) with a modestly helium-enhanced model where AY/AZ = 3. 

In Figure [5] we compare the age and metallicity derivations using the fiducial (open 
squares) and helium-enhanced (crosses) stellar models. We find that GCs which have best- 
fit metallicities higher than -1.0 dex show a lower best-fit age with the helium-enhanced 
model than the fiducial model. The use of an helium-enhanced model moves the FUV 
luminosity-weighted ages of the GCs closer to their 'true' ages, since He-rich stars are hotter 
and thus brighter in the UV than their He-poor counterparts for the same age and mass. 
Nevertheless, the modest helium-enrichment employed in the AY/AZ = 3 model does not 
lower the estimated ages sufficiently to achieve agreement with the spectroscopic ages. 



The He enrichment (AY/AZ > 90), for a sub-population of stars, adopted in 



Leeetal 



(120051 ) is significantly more extreme. To explore whether such extreme He enrichment could 



be present in the M87 GCs, we compare the FUV luminosity-weighted age of to Cen with 
its probable true age. To compute the true age of u Cen, we calculate the mass-weighted age 




Figure 6. TOP PANEL: Best-fit age and metallicity estimate for u) Cen and one-sigma contour. BOTTOM LEFT: Marginalised 
metallicity probability density function for ui Cen. BOTTOM RIGHT: Marginalised age probability density function for ui Cen. 



of the best model of 



Lee et al.l (120051 ). which accurately reproduces the observed CMD of uj 
Cen (see Table 1 in their paper). This mass- weighted age is ~ 12 Gyr. We then compute 
the FUV luminosit y-weighted age for th is GC, using its (15 — V) , (25 — V) and (V — I)o 
colours, taken from Iporman et al. (199^), and by performing the parameter estimation with 



the fiducial AY/AZ = 2 stellar model. 

As Figure [6] indicates, the best-fit metallicity for u Cen is around —1.3 ± 0.21 dex (for 
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Figure 7. Age and metallicity est imates for M87 GC s, derived using UV and optical photometry (black circles with error 
bars). Ages of M87 GCs studied bv lCohen et al.l (jl99ST ) . determined from Lick indices, are shown using the open triangles and 
grey shaded region. The FUV luminosity-weighted age (15.1 Gyr) and mass-weighted age (~12 Gyr) of ui Cen are shown using 
the large filled triangles. 



comparison Dorman et a!., 1995 quote the metallicity of u Cen as -1.6 dex) and the best-fit 
age is 15.08 ±0.4 Gyr. We find that the discrepancy between the FUV luminosity-weighted 
age and the mass-weighted age of uj Cen is ~ 3 Gyr, wh i ch is very similar to what we found 
for the M87 GC population detected in the Sohn et all J2OO6I ) study. 

We collate all our findings in Figured The filled circles show age and metallicity estimates 
for M87 GCs, derived using the AY/AZ = 2 model. Ages for M87 GCs, determined from 
the Lick indices, are overplotted using open triangles. The FUV luminosity-weighted (15.1 
Gyr) and mass- weighted (12 Gyr) ages of uj Cen are shown by the large filled triangles. The 
striking correspondence between the discrepancies in the luminosity and mass-weighted ages 
in both uj Cen and the M87 GC population, indicates that a similar He enrichment to what 
is conjectured in uj Cen would produce a similar enhancement of FUV flux in M87 GCs, 
causing their ages to be similarly overestimated by ~ 3 Gyr. 
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5 SIMULATING THE M87 GC POPULATION 



Figure [3] shows the observed UV GC data from M87 and the Milky Way. We find that 



even optica l 



Sohn et al 



y-red GCs show strong UV fluxes and blue UV — V colours (FUV — V ^ 2). 



(120061 ) found that all optically bright (V < 22) GCs in M87 are also UV bright, 



regardless of their metallicities. Both in the Milky Way and M87 samples, there appears to 
a be a positive and thus sensible correlation between V — I and FUV — V; but the M87 
GCs are brighter in the UV than their Milky Way counterparts by at least a magnitude 
or more. We note that, for a sample of M31 GCs detected in the FUV by the GALEX 
satellite mission, the same stellar models yield best-fit ages of 10-12 Gyr, using an identical 
parameter estimation technique (S. -C. Rey, private communication). 

To investigate possible reasons for these surprising features of the M87 sample, we simu- 
late a population of GCs using the GCLF and GC metallicity distribution function (MDF) 
in M87 (Figure EJ), taken from lKundu et all (jl999 ). Each GC is allocated a formation age, 
drawn from a uniform distribution of values between 9 and 13 Gyr. This yields a synthetic 
population of GCs, each described by a V^-band magnitude, a metallicity and a formation 
age. Using the fiducial stellar models (AY/AZ = 2), we then compute the apparent FUV 
magnitude of each GC, given its apparent V-band magnitude, metallicity and age. This 
results in the distribution of apparent FUV magnitudes shown in Figure [91 We now explore 
the characteristics of GCs that we wo uld expect to detec t in this simulated sample, given the 



FUV detection limit of 25 mag in the 



Sohn et al. 



(120061 ) study. Note that, although we have 



performed multiple (~50,000) simulations, the plots in this section describe results from a 
typical simulation of ~1000 objects. 

Based on a large number of simulations, we find that the detectable fraction of normal old 
GCs, is 2.1 ± 0.6 percent. We note t hat this valu e is rea sonably robust against the variation 



in the input distributions. Since the 



Sohn et al. 



(120061 ) study (excluding Field 4 because it 



has a lower detection li mit and is not r eprese ntative of the other fields) covers ~ 30 percent 



of the WFPC2 image of 



Kundu et al. 



(119991 ). we would expect it to have ~ 300 GCs in its 
field of view. In an ideal scenario, where all the detectable GCs are indeed picked up by the 
observations, a detection rate of ~2 percent translates to ~6 clusters detected i n the FUV. 



'he a ctual detection rate is significantly higher, with 66 sources det ected by 



( 120061 ). of which 50 are confirmed as GCs through comparison wit 



i the 



study. The number is lowered to 38 after cross-matching with 



Sohn et al. 



Kundu et al 



Jordan et al. 



1999] ) 



20021 ). We 
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Figu re 9. TOP: Simulated distribution of apparent GC FUV magnitudes. The solid line shows the detection limit of the 
ISohn et al.l l l2006t) study. BOTTOM: App arent FUV magni tudes vs. apparent V-band magnitudes for the simulated GCs. The 
solid line shows the detection limit of the lSohn et al.l J2006T) study. The plots show a typical simulation of ~1000 objects. 
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therefore find that the detection rate achieved by this study is up to eight times higher 
than what might be expected on the basis of canoncial stellar model predictions. Given 
the expected detection rate of ~ 2%, the observed detection rate of ~17% (50 out of 300) 
achieved by this study is a 24 sigma event! 

We finish this section by comparing the FUV and optical colours of the observed sample 
to those in the simulations (Figure [TUT) . As expected from the analysis presented in the 
preceding sections, we see a clear separation between the parameter space occupied by the 
simulated GCs and that occupie d by the observed sample - we illustrate this using grids 



which plot th e 



Kundu et al 



SohnetaL 



(|1999h and 



( 2006) FUV pho tometry against optical photometry from both 



Jordan et al. 



((2002). 



6 SUMMARY AND DISCUSSION 

We have studied the UV photometry from HST / ST IS and optical photometry from HST/WFPC2 
photometry of 38 GCs in the Virgo elliptical M87. The dataset represents the best photo- 
metric GC dataset covering both the UV and the optical ranges. Parameter estimation of 
the age and metallicity of each M87 GC, using stellar models with standard He enrichment 
(AY/AZ ~ 2), results in almost all GCs having a best-fit age of 14-18 Gyr, significantly in 
excess of the currently accepted age of the Universe. 

The metallicity estimation leads to an age distribution which shows no clear trend with 
metallicity at least in our limited sample. The marginalised metallicity PDF for the entire 
sample appears multi- modal. A KMM analysis for bimodality yields two peaks at -1.33 



Jordan et al. 


(2002 


) and 


Kundu et al. 


(1999) 



limited because our sample is a Ft/V-selected subset of the parent GC population and thus 
are affected by selection effects (and also low number statistics). The derived ages of the GCs 
in this study are significantly larger than the Lick-index-based ages of a more representative 
sample of the M87 parent GC population. 

We have presented a plausibility argument that the 'overestimate' in the derived ages may 
be driven by enhanc ed FUV f uxes fr om anomalous HB components in these GCs, similar to 
those postulated by iLee et al.l (120051 ) in the Galactic GC u Cen, where a sub-population of 
EHB stars, which are hotter and super-helium-enriched compared to the primary population 
of normal HB stars, enhance the overall FUV flux more than would be normally expected. 
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2 4 

FUV-V 




Figure 10. TOP: (FUV — V) vs. (V — I) for observed (filled circles with error bars) and simulated GCs (crosses) plotted 
on top of a model grid generated from stellar models with the fi ducial value of He enrichment (AY/AZ = 2). Filled circles 
repres ent t he 38 GCs used in this study after cross-matching the ISohn et al. I 1120061 1 cat alog with that of both 
2005 ) and lKundu et a ll jl999l). We a lso indicate, using open circles, GCs which appear in lSohn et al.l l|2006ri and 



_ , ,^99]). We a lso indicate, using open circles, GCs which appear in So hn et al . (2006) and Kundu et ah. 

19991) but not in I Jordan et alj d2002h , and which arc therefore excluded from the parameter estimation. The lowest model age 



Jordan et al 



plotted is 1 Gyr and the largest age plotted is 15 Gyr. Ages 1,5,10 and 15 Gyr are shown using filled circles of increasing sizes. 
BOTTOM: (FUV — V) vs. (F336ui — .F547m) for observed (filled circles with error bars) and simulated GCs (crosses), plotted 
on top of a model grid generated from stellar models with the fiducial value of He enrichment (AY/AZ = 2). 
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To support this hypothesis, we have first shown that replacing the fiducial stellar model, 
which assumes AY/AZ = 2, with a modestly helium-enhanced stellar model (AY/AZ = 3) 
lowers the best-fit ages of only the metal-rich M87 GCs by ~ 1 to 2 Gyr. We have then 
compared the FUV luminosity- weighted age of u Cen, derived from its UV photometry in 
(119951 ) , t o its 'true' age, c alculated from the mass- weighted mean age of the 



Dorman et al. 



best-fit u) Cen model of 



Lee et al. 



(120051 ). We have found that this "true" age is around 12 
Gyr compared to a FUV luminosity- weighted age of around 15.1 Gyr - a discrepancy of 3 
Gyr, which is very similar to the case of the FUV bright M87 GCs in our sample. 

Given the striking similarity in the discrepancy between the FUV luminosity-weighted 
and mass-weighted ages of both to Cen and M87 GCs, we suggest that a similar type of 
extreme He enrichment scenario could, in principle, produce the enhanced EHB components 
in the M87 GCs. Given the fact that our M87 sample is FUV selected, it is not surprising 
that we detect mainly those GCs which show enhanced EHB components (and therefore 
enhanced FUV fluxes), leading to overestimates of their ages when stellar models with 
fiducial values of primordial He abundance are applied. This hypothesis is supported by 
a simulation of M87 GCs, which ind icates that the de tection rate for old GCs, given the 



GCLF, GC MDF, area covered by the 



Sohn et al. 



(120061 ) study, and assuming formation ages 



between 9 and 13 Gyr, should only be around 2 percent - the actual detection rate is up 
to eight times higher. Since few GCs are expected to be observed but a significant number 
are detected, it is natural to conclude that the FUV fluxes are enhanced in the majority 
of the detected GCs. Although we cannot prove its valitidy based on this argument, the 
super-He-rich hypothesis offers a reasonable explanation. 

It is clear that drawing conclusions about the global FUV properties of the M87 GC 
population, based on this "super-helium-enriched" subsample alone, is not practical. Deeper 
observations are required to sample both more GCs with enhanced EHB components and 
also normal GCs which conform to the expectations of the stellar models for old populations 
(such as those found in M31). The presence of super-He-rich stellar populations in Galactic 
GCs such as to Cen and in the M87 GC systems raises the intriguing possibility that this 
phenomenon may be present, in varying degrees, in most systems in the nearby Universe. 



A good example of such a case could be 



dCode fc Welc 



1979 



Burstein et al. 



; he UV upturn p henomenon of elliptical galaxies 



1988 



O'Conneil 



19991 ). Indeed, this study may be the 



first to discover signatures of apparent He enrichment in extragalactic objects. Future UV 
surveys of nearby GC systems are therefore highly anticipated, both to establish the possible 
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ubiquity of the He enrichment phenomenon and to measure the degree to which stellar 
populations may be He enriched. 
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